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NUMERICAL SCHUBERT CALCULUS 


BIRKETT HUBER, FRANK SOTTILE AND BERND STURMFELS 


Abstract. We develop numerical homotopy algorithms for solving systems of 
polynomial equations arising from the classical Schubert calculus. These homo- 
topies are optimal in that generically no paths diverge. For problems defined by 
hypersurface Schubert conditions we give two algorithms based on extrinsic defor¬ 
mations of the Grassmannian: one is derived from a Grobner basis for the Pliicker 
ideal of the Grassmannian and the other from a SAGBI basis for its projective 
coordinate ring. The more general case of special Schubert conditions is solved 
by delicate intrinsic deformations, called Fieri homotopies, which first arose in the 
study of enumerative geometry over the real numbers. Gomputational results are 
presented and applications to control theory are discussed. 
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1. Introduction 


Suppose we are given linear subspaces Ki,... , of with dim Ki = m+l — ki 
and ki + ■ ■ ■ + kn = mp. Our problem is to find all p-dimensional linear subspaces of 
Qm+p meet each Ki nontrivially. When the given linear subspaces are in general 

position, the condition ki + ■ ■ ■ + kn = mp guarantees that there is a hnite number 
d = d{m,p, ki,... ,kn) of such p-planes. The classical Schubert calculus gives 
the following recipe for computing the number d. Let hi,... , be indeterminates 
with degree{hi) = i. For each integer sequence Ai > ■ ■ ■ > Ap+i we dehne the 
following polynomial: 


( 1 ) 




Here hg := 1 and hj := 0 if i < 0 or i > m. Let / be the ideal in Q[hi,... , 
generated by those Sx with m > Ai and Ap+i > 1. The quotient ring Am,p ■= 
Q[hi,... , hni]/I is the cohomology ring of the Grassmannian of p-planes in C^+p, It 
is Artinian with one-dimensional socle in degree mp. In the socle we have the relation 


( 2 ) 


d ■ hk^hk2' ■ ■ ^kn 


Thus we can compute the number d by normal form reduction modulo any Grobner 
basis for I. More efficient methods for computing in the ring Am,p are implemented 
in the Maple package SF ^ . 


In the important special case fci = • ■ ■ = = 1 there is an explicit formula for d: 

^ ^ l!2!3!---(p-2)!(p-l)!- {mp)\ 

m\ (m-|-l)! {m+ 2)! ■ ■ ■ —1)! 

The integer on the right hand side is the degree of the Grassmannian in its Pliicker 
embedding. This formula is due to 0; see also |^, XIV.7.8] and Section 2.3 below. 

The objective of this paper is to present semi-numerical algorithms for computing 
all d solution planes from the input data Ki,... , Kn- This amounts to solving certain 
systems of polynomial equations. Our algorithms are based on the paradigm of 
numerical homotopy methods 0 . 0 . 0 . 

Homotopy methods have been developed for the following classes of polynomial 
systems: 

1. complete intersections in affine or projective spaces 0, 0], 


2. complete intersections in products of projective spaces [p7| , 

3. complete intersections in toric varieties 0, 0 . 


In these cases the number of paths to be traced is optimal and equal to the standard 
combinatorial bounds: 
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1. the Bezout number (= the product of the degrees of the equations) 

2. the generalized Bezout number for multihomogeneous systems 

3. the BKK bound |^, ^ Q (= mixed volume of the Newton polytopes) 


None of these known homotopy methods is applicable to our problem, as the fol¬ 
lowing simple example shows: Take m = 3, p = 2, and fci = ■ ■ ■ = fcg = 1, that is, 
we seek the 2-planes in which meet six general 3-planes nontrivially. By formula 
(1.3) there are d = 5 solutions. Formulating this in Pliicker coordinates gives 11 ho¬ 
mogeneous equations in ten variables, the hve quadrics in display (|^) below and six 
linear equations (|^. A formulation in local coordinates (P) has 6 quadratic equations 
in 6 unknowns, giving a Bezout bound of 64. These 6 equations all have the same 
Newton polytope, which has normalized volume 17, giving a BKK bound of 17. 

In Section 2 we give two homotopy algorithms which solve our problem in the 
special case ki = ■ ■ ■ = = 1, when the number of solutions equals (1.3). The hrst 

algorithm is derived from a Grobner basis for the Pliicker ideal of a Grassmannian 
and the second from a SAGBI basis for its projective coordinate ring. (See Eg 
or 1^, Gh. 11] for an introduction to SAGBI bases). Both the Grobner homotopy 


and SAGBI homotopy are techniques for hnding linear sections of Grassmannians in 
their Pliicker embedding. 

In Section 3 we address the general case of our problem, that is, we describe a 
numerical method for solving the polynomial equations dehned by special Schubert 
conditions. This is accomplished by applying a sequence of delicate intrinsic deforma¬ 
tions, called Fieri homotopies, which were introduced in |^. Pieri homotopies hrst 
arose in the study of enumerative geometry over the real numbers |^. For the 
experts we remark that it is an open problem to hnd Littlewood-Richards on homo¬ 
topies, which would be relevant for solving polynomial equations dehned by general 
Schubert conditions. 

A main challenge in designing homotopies for the Schubert calculus is that one 
is not dealing with complete intersections: there are generally more equations than 
variables. In Section 4 we discuss some of the numerical issues arising from this 
challenge, and how we propose to resolve them. In Section 5 we discuss applications 
of these algorithms to control theory and real enumerative geometry. Finally, in 
Section 6 we present computational results. 

In closing the introduction let us emphasize that all homotopies described in this 
paper are optimal in the sense the number of paths to be traced equals the number 


d. This means that for generic input data Ki 


Kn no paths diverge. 


2. Linear equations in Plucker coordinates 

The set of p-planes in Grass{p,m + p), is called the Grassmannian of p- 

planes in This complex manifold of dimension mp is naturally a subvariety of 

the complex projective space Pv p / .To see this, represent a p-plane in C^+p as 
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the column space of an [m + p) x p-matrix X = {xij). The Plucker coordinates of 
that p-plane are the maximal minors of X, indexed by the set of sequences 

a ■. 1 < a\ < a 2 <■■■< Oip < m + p\ 


( 4 ) 



- 1 


Xqi p 


^OLpP . 


This section deals with the = 1” case of the problem stated in the Introduction. 
Given mp general m-planes Ki,... , Kmp, we wish to hnd all p-planes X which meet 
Ki,... ,Kmp nontrivially. This geometric condition translates into linear equations 
in the Pliicker coordinates: Represent X as an [m + p) x p-matrix as above, represent 
Ki as an (m + p) x m-matrix, and form the (m + p) x (m + p)-matrix [X \ Ki]. Then 


X n Ki ^ {0} if and only if det [X | iPj] = 0. 


Laplace expansion with respect to the hrst p columns gives 


( 5 ) 


det[X|X,] = Ci-[a], 

„e([-+Pl) 


where is the correctly signed maximal minor of complementary to a. Hence 
our problem is to solve mp linear equations (|^) on the Grassmannian. The number of 
solutions is the degree of the Grassmannian in its Plucker embedding, which is given 
in (1.3). 

The Grassmannian is represented either implicitly, as the zero set of polynomials 
in the Plucker coordinates, or parametrically, as the image of the polynomial map 
(§). These two representations lead to two different numerical homotopies. The im¬ 
plicit representation gives the Grobner homotopy in Section p.2| while the parametric 


representation gives the SAGBI homotopy in Section |2.3| . The first is conceptually 
simpler but the second is more efficient. In both methods the number of paths to be 
traced equals the optimal number in (H). 


2.1. An example. We describe the two approaches for the case (m,p) = (3, 2). The 
Grassmannian of 2-planes in C® has dimension 6 and is embedded into P®. Its degree 
(|^) is hve. The Grobner homotopy works directly in the ten Plucker coordinates: 

[12], [13], [14], [15], [23], [24], [25], [34], [35], [45]. 
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The ideal of the Grassmannian in the Pliicker embedding is generated by hve 
quadrics: 


( 6 ) 


[14]|23] - |13][24] + [12][34], 

|15]|23] - |13][25] + [12][35], 

|15||24| - |14]125] + 112|[45|, 

|15||34| - |14]135] + [1311451, 

[25|[34| - [24][35] + [231[451. 


This set is the reduced Grobner basis for Js 2 with respect to any term order which 
selects the underlined terms as leading terms (see Proposition 2.1 below). 

Our problem is to compute all 2-planes which meet six sufficiently general 3-planes 
7^1,... , Kq nontrivially. This amounts to solving (P) together with six linear equa¬ 
tions 


,, . [12| + q3. [13] + q^. [14] + q^ ■ [is] + ■ [23] 

+ G 24 ■ [24] -|- (725 ■ [25] + (^34 ■ [34] + (Pgg • [35] -|- ■ [45] = 0, 

for i = 1,... , 6. This is an overdeter mined system of 11 equations in 10 homogeneous 
variables. To solve it we introduce a parameter t into (P) as follows: 


[14] [23] - t ■ [13] [24] + . [12] [34] 

[15] [23] - f2.[i3][25] + t4.[i2][35] 

[15] [24] - t ■ [14] [25] + ■ [12] [45] 

[15] [34] - . [14] [35] + ■ [13] [45] 

[25][34] - t -[24][35] + f2.[23][45] 


0 , 

0 , 

0, d) 

0 , 

0 . 


We call (P') the Grobner homotopy because this is an instance of the flat deformation 
which exists for any Grobner basis; see |^, Theorem 15.17]. The flatness of this 
deformation ensures that, for almost every complex number t, the combined system 
(d)&(0) has hve roots. 

For t = 0 the equations (^) are square-free monomials. We decompose their ideal: 

( 8 ) 

{ [141[231, [1511231, [1511241, [15][341, [251[34] > 

= (123], |24], |34]) n {[15], [23], [34] > n ([15], [23], [25]) 

n {[14], [15], [34] > n ([14], [15], [25] >. 


The hve distinct solutions for f = 0 are computed by setting each listed triple of 
variables to zero and then solving the six linear equations (0) in the remaining seven 
variables. Thereafter we trace the hve solutions from f = 0 to f = 1 by numerical 
path continuation. At f = 1 we get the hve solutions to our original problem. 
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We next describe the SAGBI homotopy. For this we choose the local coordinates 


1 

0 

2:21 

2:22 

2:31 

2:32 

2:41 

2:42 

0 

1 


on the Grassmannian. Substituting X into (0) we get six polynomials in six un¬ 
knowns: 

fq'i Qa ■ (2^212:32 - 2:222:31) + C'k ■ (2:212:42 - 2:222:41) + <^34 ■ (2^312:42 - 2:322:41) 

+ ^25 ■ ^21 + C35 • 2:31 -|- G45 • 0:41 -|- CI2 ■ 2:22 + GI3 • Xs 2 + GI4 • 0:42 + C\^. 

To solve these six equations we introduce a parameter t as follows: 

<^23(^212:32 - t ■ X22X31) + ^24(0:212:42 - ■ 2:222:41) + G34 (0:310:42 - t ■ X32X41) 

+^25 ■ 2:21 + G35 • 0:31 -|- G45 • 0:41 -|- CI2 ■ X22 + GI3 ■ 0:32 -|- GI4 • 0:42 -|- GJ5. 


The system (^') has five complex roots for almost all t G C. For t = 0 we get a 
generic unmixed sparse system (in the sense of |^) with support 


A — {1, X 21 , X 22 , 2 : 31 , 0 : 32 , 0 : 41 , 0:42 , 2 : 210:32 , 0 : 210:42 , 0 : 310:42 }. 


We identify A with a set of ten points in Z®. Their convex hull conv{A) is a 6 - 
dimensional polytope with normalized volume five. We can therefore solve (§') for 
t = 0 using the homotopy method in or |^, provided the input data Ki,... ,Kq 


are sufficiently generic. Tracing the hve roots from f = 0 to f = 1 by numerical path 
continuation, we obtain the hve solutions to our original problem. 


2.2. Grobner homotopy. We next describe a quadratic Grobner basis for the dehn- 
ing ideal of Grass{p,m + p). Let S be the polynomial ring over C in the variables [a] 
where a G We dehne a partial order on these variables as follows: [a\ < [j3] 

if and only if a, < f3i for i = 1,... ,p. This partially ordered set is called Young’s 
poset. Figure |] shows Young’s poset for {m,p) = (3,2). 

Fix any linear ordering on the variables in S which rehnes the ordering in Young’s 
poset, and let -< denote the induced degree reverse lexicographic term order on S. 
Let Im,p be the ideal of polynomials in S which vanish on the Grassmannian, that 
is, /m,p is the ideal of algebraic relations among the maximal minors of a generic 
{m + p) X p-matrix X. The Grobner homotopy is based on the following well-known 
result. 

Proposition 2.1. The initial ideal in^{Im,p) is generated by all guadratic monomi¬ 
als [a][/5] where ai < j3i and aj > /3j for some i,j G {1,... ,p}. 
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[ 12 ] 

[ 13 ] 

/ \ 

[ 14 ] [ 23 ] 

/ \ / 

[ 15 ^ 

[ 25 ] [ 34 ] 

\ / 

[ 35 ] 

[ 45 ] 

Figure 1. Young’s poset for (m,p) = (3,2). 

In other words, is generated by products of incomparable pairs in Young’s 

poset. Let Cm,p denote the set of all maximal chains in Young’s poset. For example, 

C3,2 = { {[12], [13], [14], [15], [25], [35], [45]} , {[12], [13], [14], [24], [25], [35], [45]} , 

{[12], [13], [14], [24], [34], [35], [45]} , {[12], [13], [23], [24], [34], [35], [45]} , 

{[12], [13], [23], [24], [25], [35], [45]}} 

A standard result in combinatorics [^ states that the cardinality of Cm,p equals the 
number (j^). From Proposition we read off the following prime decomposition 
which generalizes (§): 

(10) in^{Im,p) = n ^ 

CGCm,p 

For a proof of Proposition see |]^ §XIV.9] or ||, Theorem (4.3)] or 
In these references one hnds an explicit minimal Grobner basis for Im,p, which is 
classically called the set of straightening syzygies. In the special case p = 2 the 
straightening syzygies coincide with the reduced Grobner basis: 

Proposition 2.2. If p = 2 then the reduced Grobner basis of Im,p consists of the 
three-term Pliicker relations [7/][fcj] — [7fc][j7] + [7j][fc/] where l<i<j<k<l< 
m + p = m + 2. 

For p > 3 the straightening syzygies and the reduced Grobner basis do not coin¬ 
cide, and they are complicated to describe. For our purposes the following coarse 
description suffices. Let Std be the set of all quadratic monomials in S which do not 
he in in^{Im,p)- The reduced Grobner basis consists of elements of the form 

(11) HI/31 - E Y’LWM 

[ 7 ] [(5]eStd 

where [a][/3] runs over all generators of m^(/m,p). The constants are integers 

which can be computed by substituting (H) into ([Tl| ) and solving linear equations. 


38, §3.1]. 
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The term order -< can be realized for the ideal Im,p by the following choices of 
weights. We define the weight of the variable [a] = [aia 2 ■ ■ ■ ctp] to be 

( 12 ) := ^ {aj-ai-lf. 

If we replace each variable [a] in (|Tl|) by and clear f-denominators afterwards, 

then we get the Grobner homotopy: 

- E (0) 

[7][(5]g5i(i 


It can be checked that all exponents + vs — Va — vg appearing here are positive 
integers. The special case {m,p) = (3,2) is presented in (|^). 

In the Grobner homotopy algorithm, we hrst solve systems of mp linear eqnations, 
one for each chain C G Cm,p- These systems consist of the mp eqnations (^), one for 
each Ki, and the — mp — 1 eqnations 

[a] for [a] ^ G, 


snggested by the prime decomposition (|I0D . Once this is accomplished, we trace each 
of these d solntions from f = 0 to f = 1 in the Grobner homotopy (^). 

Glearly, the weights Va of (|T^ are not best possible for any specific valne of m and 
p. Smaller weights can be fonnd nsing Linear Programming, as explained e.g. in the 
proof of Proposition 1.11]. Another method wonld be to adapt the “dynamic” 
approach in to onr sitnation. This is possible since the Grobner basis in (pd]) 


IS 


reverse lexicographic: hrst deform the lowest variable to zero, then deform the second 
lowest variable to zero, then the third lowest variable, and so on. 


2.3. SAGBI homotopy. Let X = {xij) be an (m+p) xp-matrix of indeterminates. 
We identify the coordinate ring of the Grassmannian with the C-algebra generated 
by the p x p-minors of X. Gall this algebra R and write [a] {xij) for the minor indexed 
by a. Reinterpreting classical resnlts in [^, §XIV.9], it was shown in [^, Theorem 
3.2.9] that these generators form a SAGBI basis with respect to the lexicographic 
term order indnced from xn > Xi 2 > ■ ■ ■ > Xip > X 21 > ■ ■ ■ > Xm+p,p- This means 
that the initial algebra C[m>(/) : f E R] is generated by the main diagonal terms 
of the p X p-minors, 

(13) Qcr] (Xjj)) ^a 2,2 ^ 03,3 . 

The resnlting flat deformation can be realized by replacing Xij with for 

f 0 in the matrix X. If we expand [a]{xijt^'^~^'^^^~^'^) as a polynomial in t, then 
the lowest term equals times the main diagonal monomial (p!3[) , where 

p 

Wa ■= ^{aj-l){p-3). 

i=i 
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In what follows we multiply that polynomial by t 
algebra 


Rt 


C 




For any t G C consider the 

« e (t’-yl) ■. 


Then Ri is the coordinate ring of the Grassmannian, and Rq is the algebra generated 
by the monomials (|^. This is a flat deformation of C-algebras; see |]^ and |^, 

§ii|. 

The SAGBI homotopy is the following system of mp equations: 


(14) ^ G; ■ = 0 (* = 1,... ,mp). 


We reduce the number of variables to mp by introducing local coordinates as follows: 
Xu = 1 for i = 1,... ,p and Xij = 0 for i < j or i > m + j. Our original problem is 
to solve the system (|^) for t = 1. 

The flatness of the family of algebras Rt guarantees that the system (0) has the 
same hnite number of complex solutions (counting multiplicities) for almost every 
t E C. For f = 0 we get a system of linear equations in Rq: 


(15) 


• Xat,lXa2,2Xa3,3 ' ' ' Xap,p (* = !,..., mp). 

„e([-+Pl) 


In order to solve these equations we apply the symbolic-numeric algorithm in ||17|| , 
while taking advantage of the following combinatorial structures described in ||39| , 
Remark 11.11], The common Newton polytope of the equations (1^) equals the order 
polytope of the product of an m-chain and a p-chain (Sturmfels, Remark 11.11). We 
have the following combinatorial result. 


Proposition 2.3. The following five numbers coincide: 

• the right hand side of (Qj, 

• the number of linear extensions of the product of a m-chain and a p-chain, 

• the number of maximal chains in Young’s poset, 

• the normalized volume of the order polytope, and 

• the number of roots in (C*)™^’ of a generic system 


21 


The equality of (4) and (5) is a special case of Kouchnirenko’s Theorem 
as all the equations have the same Newton polytope. The order polytope has a 
distinguished unimodular regular triangulation with simplices indexed by the chains 
in Young’s poset. This regular triangulation is induced by the system of weights 
given in (|^). We may use these weights to dehne a numerical homotopy for hnding 
all isolated solutions of (p!5D. Once this is accomplished, we trace these roots from 
f = 0 to f = 1 in the homotopy (|^). 
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3. Special Schubert conditions 

We first describe a purely combinatorial method for computing the number d of 
solution planes. Instead of the algebraic relation (1.2) we shall make use of Young’s 
poset which was introduced in Section 2.2. A cover [a] < [[3] in Young’s poset deter¬ 
mines a unique index j = j{a, j3) for which 

aj + I = (dj and a, = jdi for i ^ j. 

A chain [a°] < [a^] < ■ ■ ■ < [a^] in Young’s poset is increasing at i if either f = 1, 
or else i> 1 and For instance, [12] < [13] < [14] < [24] is 

increasing at 1 and 2, but decreasing at 3. 

Given positive integers tq, ... , r^, the Fieri tree T(cq, ... , r^) consists of all chains 
of length ro + ■ ■ ■ + La in Young’s poset which start at the bottom element [1,2,... , p] 
and which increase everywhere, except possibly at lq -|- 1, Tq -|- ri -|- 1,... , tq + ■ ■ ■ + 
ra-i + 1. Here we include all initial segments of such chains and we order the chains 
by inclusion. Label a node in the Fieri tree by the endpoint of the chain which that 
node represents. Then the sequence of labels from the root to that node is the chain 
which that node represents. For example, here is T(2, 2) when m = 5, p = 2: 

[34] [25] [16] 

\ / / 

[24] [15] 

\ / 

[14] 

[13] 

[12] 

To compute the number d, partition the integer sequence ki,... ,kn into three 
parts To, ... , Tq, Tq, ... , r'^,, and q. Then d is the number of pairs {R, S) where R is 
a leaf of T (lq, ... , Tq), S is a leaf of T(rQ,... , r'^,), and the endpoints [a] of R and 
[a'] of S satisfy Pieri’s condition-. 

(16) a(<m-|-p-|-l — ap< 02 <■■•<«(, <m-|-p-|-l — 

Call this set of pairs Sols = Sols{ro ,... , r^; rg,... , r(,,). 

For instance, d = 6 for the sequence 2, 2, 2, 2, 2 with {m,p) = (5, 2): Of the 9 pairs 
{R, S) of leaves of T(2, 2), only the following 6 satisfy (pTl): (Here we represent a leaf 
by its label.) 

(17) {([25], [34]), (|34],|25]), (|25|,|25|), (|25|,|16|), (116|,125|), ([16], [16])}. 

This combinatorial rule for the number d gives the same answer as the algebraic 
rule (1.2) because the Fieri tree and Fieri’s condition (pT]) encode the structure of the 
cohomology ring Am,p with respect to its Schur basis {S\ | p > Ai > • ■ ■ > Ap+i — 0}; 
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see 


p] or 1^, §9.4], Specifically, 

n 


,ra) ■ 

^—0 AiH -hAp=roH- \~'f'a 

where i^A,(ro,...,ra) is the number of leaves in T(ro,... ,ra) with label 

[a(A)] := [Ap + 1, Ap_i + 2 ,... , Ai + p]. 

The numbers Kx,(ro,...,ra) called Kostka numbers. In Am,p we calculate 


n'*‘. = 'll'*-- 




i=l 


,2 = 0 


0=0 


-^A,(ro,... T^/) ' ' Sp ■ kg. 

X,li 

We evaluate this expression with Pieri’s formula (Proposition |3.1| below): If Ai + - ■ ■ + 
Ap + /ii + ■ ■ ■ + /ip + g = mp, then Sx ■ Sp ■ hg is either {hmY or 0 depending upon 
whether or not 

\p < m — Pi < X 2 < ■ ■ ■ < Xi < m — pp. (PTf) 

The condition ([l6f ) is equivalent to (IT) under the transformation A <->■ [a(A)]. 

These methods correctly enumerate the p-planes which meet each Ki,... , non- 
trivially because, under the isomorphism between Am,p and the cohomology ring of 
Grass{p,m + p), the indeterminate hfc. corresponds to the cohomology class Poincare 
dual to klxi, the set of p-planes which meet Ki nontrivially. Moreover, {hmY repre¬ 
sents the class dual to a point. The Pieri tree models certain intrinsic deformations 
(described in §3.3 and § p.5|) of the Grassmannian which establish this isomorphism, 
and which we shall use for computing the p-planes which meet each iPi,... , 
nontrivially. 


3.1. Basics on Schubert varieties. For vectors /i,... , fj in Q'^+p^ let (/i,... , fj) 
be their linear span. Fix the columns Ci,... , e^+p of the identity matrix as a standard 
basis for Q^+p. For a G , dehne G by aj := m -|- p -|- 1 — ap+i_j. 

A sequence a G determines a Schubert variety 

kla '■= {K G Grass{p, m + p) \ dimX fl (ci,... , e^j) > j for 1 < j < p}. 

This variety has complex codimension |a| := Oi — 1 -|- 02 — 2 -|- ■ ■ ■ -f Op — p. Similarly 
dehne 

G'q, := {X G Grass{p,m+p) \ dimXn(ea^.,... , e^+p) > p+l-j for 1 < j < p}. 

For a linear subspace N of C^+p of dimension m -|- 1 — g, dehne the special Schubert 
variety 

Gat := {X G G'rass(p, m-I- p) | dimXnX>l}. 



12 


BIRKETT HUBER, FRANK SOTTILE AND BERND STURMFELS 


This has codimension q. U N = {ei,, Cm+i-q), then Qn = ^[i,...,p-i,p+q]- 

The special Schubert variety is cut out by the system of polynomial 

equations: 

(18) X G hi AT all maximal minors of [X | iV] are zero, 


where X G Grass{p, m + p) is represented by a (m + p) x p-matrix. The Laplace 
expansion of these equations in terms of the Pliicker coordinates of X define hi at as a 
subscheme of Grass{p,m + p). These equations are redundant: select m+1—g rows of 
[X I X] such that the corresponding maximal minor of N is invertible. Consider the 
set of maximal minors of [X | X] which cover all the rows selected. This gives 
polynomial equations which generate the same ideal as all minors of [X | X]. 

For a purely set-theoretic (but not scheme-theoretic) representation of Qjsf a further 
substantial reduction in the number of equations is possible using the results of |^. 

An intersection Y n Z of subvarieties is generically transverse if every component 
of Y n Z has an open subset along which Y and Z meet transversally. In this case 
the following identity in the cohomology ring holds: 


[y n z| = [y| • [z|, 

where [W] denotes the cycle class of a subvariety W. By Kleiman’s Transversality 


Theorem [|^, subvarieties of Grass {p,m + p) in general position meet generically 


transversally. Transversality and generic transversality coincide when Y CiZ is hnite. 


Proposition 3.1 (Hodge and Pedoe, 1952, Theorem III in §XIV.4). 

Let a, a' G with |Q;|-|-|Q;'|-|-g = mp and let N be a linear sub space with 

dimension m + \ — q none of whose Pliicker coordinates vanish. Then the intersection 

(19) fla n Vl'^, n ffjv 

either is transverse consisting of a single p-plane or is empty, depending upon whether 
or not & holds. 

Proof and Algorithm: The intersection C is nonempty if and only 
if a' < a'j for j = 1,... ,p. These are the weak inequalities in (|^). We shall 
assume that they hold in what follows. The p-planes in Hq, Pi are represented by 
(m -|-p) X p-matrices X = (t^) such that 

(20) Xjj = 0 for i < a'j or a'^ < i. 

Consider the nonzero coordinate subspaces Cj := {Ca'.,... ,ea^), set C := Ci + 

■ ■ ■ I- Cp, and note that p -|- g = — dim{C). From (^) we see that 

X G Hq, n implies X C C and hence X ft X C X fl C. Therefore the triple 
intersection ([I^) is nonempty only if the following equivalent conditions hold: 

dim{Cr]N)>l dim{C)=p + q the sum C = Ci -I — -YCp is direct 
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for j = 1,... ,p — 1 (3.1) holds. 

In this case we determine (7 fl by computing vectors Qj G Cj such that C fl iV = 
{di ® 92 ® ® dp)■ (This computation is the “algorithm” part in this proof.) The 

desired p-plane X satisfies C* fl iV = X (1 N, and, in view of (p0|), this implies 
X = {gi,... ,gp). Transversality of (|I^ is verified in local coordinates for fl VL'^, 
by considering p + g — 1 independent linear forms which vanish on N. □ 

For a G define A(q!) by X{a)j := aj —j for 1 < j < p. Then Sx(a) represents 

the cycle class of (equivalently, of !!(,). If dim N = m+ I — q, then hq is the cycle 
class of Qn- Suppose |a| + la'I + q = mp. Then Proposition |3.1| implies the following 
identity in Am,p- 


*S*A(q:) ■ Sx[a') ■ hq 


{hmY if (H) holds 

0 otherwise 


This identity implies (via Poincare duality) that 


Sx(a)-hq - y^^Sx(f3), 

the sum over all (3 with |/d| = |a| + g for which satisfy Pieri’s condition (pTl). 

Call this set a * q, which is also the set of endpoints of increasing chains of length g 
in Young’s poset that begin at a. 

This last form has geometric content. In [^, explicit deformations were given that 
transform the irreducible intersection hlo, fl hlAr into the cycle Moreover, 

the branching of the components of the cycles in these deformations reflects the 
branching among these increasing chains above a. This process may be iterated to 
transform an intersection of several special Schubert varieties into a sum of triple 
intersections of the form (^), indexed by pairs {R, S) G Sols. From this sum, we 
obtain a set of start solutions indexed by Sols. Also, every intermediate cycle in 
these deformations consists of the same number (counting multiplicities) of p-planes. 
The Fieri homotopy begins with one of the start solutions and uses numerical path 
continuation to trace the sequence of curves defined by these deformations which 
connect that start solution to a solution of the original problem. 


3.2. Fieri homotopy algorithm. Given linear subspaces Ki,... ,Kn in general 
position with dim Ki = m + 1 — R and ki + ... + kn = mp, first partition Ki,... , Kn 
into three lists: 

To, • • • , La-, Tq, ... , N 

where dimLj = m + 1 — r*, dimL' = m + 1 — r', and dimA^ = m + 1 — g. 
Construct the Fieri trees T(ro,... ,ra) and T(rg,... and form the set Sols. 

Change coordinates so that Lq = (ei,... , Cm+i-ro) and L'q = ... , Om+p)- Set 

r := maxjri H-h r^, -h r'^,}. 

Given a chain R in the Fieri tree and a positive integer k, let R{k) be the fcth 
element in that chain, or, if k exceeds the length of R, then let R{k) be the endpoint 
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of R. For each (i?, S) G Sols and k from r to 0 we shall construct (in Dehnition |3^ 
below) one-parameter families and of pure-dimensional subvarieties of 

Grass{p,m + p) with the following properties: 


1- ZR,kit) c Qfi(^ro+k) and Zgf.(t) C 

2. For t = 0 or 1 and each k, Zn^kit) fl Z'^^it) fl is transverse and 0-dimensional. 
3- Zn^T-^t) = D.jK^ro+T) and Zg,j.{t) = • 

4. ZA;,fc+i(l) is a component of Z^^iQi). Likewise, is a component of 

Z's0)- 

5. Za;’o( 1) = n ■ ■ ■ n and n ■ ■ ■ n 

Property 4 is a consequence of Proposition |3.4| , the others follow from the as¬ 
sumption of genericity and the dehnition (Dehnition |3.2|) of the families Z^^kit) and 

By 2, the family W(Fi^s),k(t) over C whose hbre at general t (including f = 0 and 
t = 1) is 

W{R,s),k{t) ■= ^R,fc(^) n Zg j^{t) n Qn 

consists of a hnite number of curves. In fact, for general t (including f = 0 and f = 1), 
W(R,s),k{t) has the following general form (see Dehnition for the precise form): 


^{R,s),k{t) — Dq, n n Dmi n • ■ ■ n Dms, 


where Mi,... , are linear subspaces with Mg = N, which depend upon R, S, k, 
the subspaces Li,... , La, L[, ... and at most two of the Mj depend upon t. 
Also, a and a' depend upon R, S, and k with the typical case being a = R{ro + k) 
and a' = S{r'Q + k). 

The numerical homotopy dehned by the curves W(^R^s),kif) may be expressed in a 
parameterization X = {xi^j) of an open subset of fl 

(21) Xij = 0 if i < a'- or aj < i and x^.j = 1, 

where 5 := S{r'Q + t). The equations for W(^R^s),kit) me then 

maximal minors [X | Mj] = 0 i = 1,... , s. 


The curves of W(R^s),k(t) dehne the sequences of homotopies in the Pieri homotopy 
algorithm as follows: For {R, S) G Sols, let X(^r^s),t be the (unique by Proposition p.l|) 
p-plane in flDjv = W^r^s)A^)- % 3 and 4, X(^r^s),t e IF(ii,s),r-i(0) 

and hence lies on a unique curve in W(R^s),T-iit)- Use numerical path continuation 
to trace this curve from f = 0 to f = 1 to obtain which is a point of 

fU(ij,s),r_2(0), by 4. Then X(^r^s),t-i lies on a unique curve in IF(R,s),r- 2 (t), which 
we trace to hnd X(^r^s),t -2 G 2 ( 1 )- Continuing this process, after tracing r 

curves, we obtain ^ hF(R, 5 ),o(l), which is a solution to the original system, 

that {X(K,s),o I {R, S) G Sols} consists of 

all the solutions to the original system. 


by 5. We show shall prove in Theorem 3.5 
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3.3. Definition of the moving cycles The cycle will depend upon 

the choice of a general upper triangular {m + p) x (m + p)-matrix F with I’s on its 
anti-diagonal, 


the fcth link in the chain R, and the data Li,... ,La. The key ingredient of this 
dehnition of Zn^kit) is the construction of a one-parameter family of linear subspaces 
Aj(t) in Dehnition |3.3| , which depends upon F. The matrix F is hxed throughout 
the algorithm, its purpose is that {em+p-j, ■ ■ ■ ,em+p) equals the span of the hrst j 
columns of F, and these columns are in general position with ei,... , Cm+p- The subtle 
linear degeneracies of Aj(t) as t 0 are at the heart of this homotopy algorithm, as 


well as the explicit proof of Pieri’s formula Theorem 3.6], which we state below 
(Proposition [3^). 


Definition 3.2. 


1. If ri H- \-ra< k, then set Za^kit) := Qai^ro+r)- 

2. Otherwise, dehne c by ri ■ -f Vc-i < k < ri + ■ ■ ■ + Vc, and set i ■= k — ri — 

-Lc-i, a := i?(ro -h ri H-h rc_i), and (3 := i?(ro -F k). 

(a) If i > 0 and (3p > Up, then /3 -|- (0,... , 0, Tc — i) = R{rQ rc), and we 

set 


^R,k{t) •= n ^Lc+1 n ■ ■ ■ n ki^a- 

(b) Otherwise, let Aj(t) be the 1-parameter family of linear subspaces given by 
Dehnition B.3|, where we let L := Lc and r := Tc- 


If i = 0, then Ao(l) = Ac, a = (3 = R{ro + k), and we set 


ZR,k{t) ■= Dij(ro+fc) n DaoP) 0 n • ■ ■ n 

If i > 0, let j be maximal such that j3j > aj. Then j < p as j = p is case 
2(a). Set 


^R,k{t') . ^ji(^ro+k) 0 ,egv ) 0 O ■ ■ ■ 0 

We dehne Z'gj^{t) similarly, but with the matrix F replaced by a lower triangular 
matrix with I’s on its diagonal, and (ei,... , replaced by ... , em+p)- 


Definition 3.3. Let F be an upper triangular matrix with I’s on the anti-diagonal 
and L be a general (m -|- 1 — r)-plane, represented as a (m -|-p) x (m -1- 1 — r)-matrix 
with columns /i,... , Im+i-r- Construct a (m -|- p) x p-matrix U = (ui,... , Up) as 
follows: Reverse the last m + p — ap columns of F, then remove the columns indexed 
by a]^,... ,a^. 

For each 0 < i < r, dehne a one-parameter family of (m + p) x (m-1-1 — r)-matrices 
Ai(t) for t G C as follows: 
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1. If z = 0, then the bth column of Ao(t) is f ■ 4 + (1 — f) ■ Ub- 

2. For 0 < z < r, the 6 th column of Aj(f) is 

t ■ Ub+i-i + {1 -t) ■ Ub+i b + i- l<ap-p 

^ 1 (f ^) " ^p+l-\-i—r 6 -f* z 1 cyp p 

Ub+i-1 b + i — l>ap — p 

3.4. An example. We give an example illustrating these definitions and the Fieri 
homotopy algorithm. Let Lq, Li, Lq, L[, and N be general 4-planes in C^. We give a 
sequence of homotopies W(^ji^s),k(t) for fc = 2 , 1,0 for finding one of the six 2 -planes 
which meet each of the five given 4-planes nontrivially. 

Here, {m,p) = (5, 2 ) and ki = ■■■ = = 2 so that r = 2 . Construct the set Sols 

as in (Q). Let {R, S) G Sols be the following two sequences: 

R := [ 12 ] < [13] < [14] < [24] < [25], S := [ 12 ] < [13] < [14] < [15] < [16]. 

Let Cl,... , ey be the columns of a 7 x 7-identity matrix, a basis for C^. Suppose that 
Lq = (e 4 , 65 , 66 , 67 ) and Lq = (ci, 62 , 63 , 64 ) and represent Li and L[ as 7 x 4-matrices. 
Then hi[ 14 ] fl fl fl H Qn is the set of 2-planes which meet all five linear 
subspaces nontrivially. 

We first hnd the plane ^ ^[ 25 ] nhlAr, using the algorithm in the proof 

of Proposition ^T|. Suppose that N has the form 

n 

~T\ ’ 

where I is the 4x4 identity matrix, and rz is a 3 x 4-matrix. In this case, Ci = 
( 61 , 62 , 63 ) and C 2 = ( 60 ), hence C = ( 61 , 62 , 63 , 66 ). Thus the intersection C fl A is 
generated by the third column of A, and so 5 ) 2 is represented by the matrix: 




Following Definition p. 3 |, we have: 


■ rzi3 O' 
^^23 0 
^^33 0 

0 0 
0 0 
0 1 
0 0 


Ao(t) = {tli +{l-t)uiRl2 +{l-t)u2Rh +{l-t)uzRh + 

Ai(t) = (tzzi - 1 - (1 - t)zZ2,tZZ2 + (1 - t)M5,M3,ZZ4). 


A((t) is defined similarly. 
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We describe the families W(^R^s),kit) for /c = 2,1, 0 in local coordinates for 
determined by the sequence [16]: 


1 

0 ■ 

X21 

0 

X31 

0 

X41 

3^42 

0 

X52 

0 

1 

0 

X72 . 


The family W(^ji^s), 2 it) is the constant family 5 ) 2 } = 1125 hi fl Assuming 
that 77,13, which is the [1457]th Pliicker coordinate of N, is non-zero, then X(^r^s ),2 
may be expressed in these local coordinates: 


X{R,S ),2 


1 0 ■ 
7123/7113 0 

nss/nis 0 
0 0 
0 0 
0 1 
0 0 


When k = 1, hrst consider the dehnition of ZR^i{t). Here we are in case 2(b) 
with 13 = [24] and z > 0, so that 13'^ = [46]. Since Ai(f) C (ei,... ,60), we have 
ZR,i{t) = 11[24] n HaiP). For the definition of we are in case 2(a), so that 

Z^ i(f) = Hjig]. Hence 

W(R,s},i(t) = 11[24] n llAi(t) n Hjie] hi Hat. 

This has 3 linear equations 0:42 = X 52 = X 72 = 0, which describe fl[ 24 ] fl Hjig], and 7 
non-trivial equations, the vanishing of the maximal minors of [X \ Ai(f)] and [X | A^], 
which describe HaiP) H Hat. 

For fc = 0, 7 = 0 and we are in case 2(b) for both ZR^it) and Z'^^it) so that 

W(R,s),o{i) = ii[i4] n HaoP) n Iiji4] n iiA^(t) n Hat- 

This has 21 non-trivial equations, the vanishing of the maximal minors of [X \ Ao(t)], 
[X I A(,(f)], and [X I X]. 


3.5. Proof of correctness. We describe the Fieri deformations linking the families 
Zr, k{t) for k from ri-f- ■ --l-rc-i to ri rc, which establishes Property 4 of 

. We also show that the set {X(a;,s),o I {R, S) G Sols} consists of all 
the solutions to the original system. 


ZR,k{t) 


m 


53.2 
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Consider the dynamic part of Zr, k{t), namely whichever of 

^R(roA -h^c)) ^i?(ro+fc) C hlAo(i)) ^-R(ro+fc) C hlAi(pn(ei,eov }) 

^p+i-j 

appeared in the dehnition of We call this cycle where L ■= and 

a = R{ro + ■ ■ ■ + Fc-i), and /3 = Rir^ + k). 

For j3 G a*i and 7 G a * (i + 1) write /3 -<q, 7 if 7 covers /3 and j{/3, 7 ) > j{a, (3) : = 
max{j I Pj > aj}. This partitions a * (i + 1) into sets {7 | /3 7 } for /3 G a * i. 

Proposition 3.4. [^, Theorem 3.6] 

Let a, L, Ai{t), andYa^js^Lit) he as above. Then 

1. For all t, hlo, n hlAo(t) is generically transverse. 

2 . Ya 3 Rt) is free of multiplicities for all t and irreducible for t ^ 0. 

3. Ifi’^r-1, then Pa,/ 3 ,L( 0 ) = ■ 

4. If /3 e a * (r - 1), then Ta,/ 3 ,L( 0 ) = E/ 3 -<, 7 ^ 7 - 

By 3, the cycle class of is independent of i and t, and it equals the 

cycle class of fl hlAo(i) = H ^2^. By 4, we see that the cycle classes of hla fl 12^ 
and E/ 3 eo*r coincide, furnishing another proof of Pieri’s formula. Property 4 of 
Zr, k{t) follows from assertion 3. 

Theorem 3.5. When Ki, ... ,Kn are generic, the Fieri homotopy algorithm finds 
all p-planes which meet each Ki, ... , Kn nontrivially. That is, 

{E-R:'S'),o I {R, s) G Sols} = Qki n ■ ■ ■ n 

Proof. Note that for any i?, S' G Sols, the families ZR^k{t),Z's^f.{t), and fF(R, 5 ).fc(^) 
depend only upon the initial segments i?(0),... , R{ro + k) and S'(O),... ,S{r'Q + k) 
of R and S. 

By construction, the original system fl ■ ■ ■ fl VLki coincides with hF(_R, 5 ),o(l )5 
for any (i?, S) G Sols. We inductively construct chains R G T(ro,... ,ra) and S G 
T^Tq, ... , r^/), and p-planes for 0 < k < r such that 

Xk+i G hF(R^5)_fc(0) n hF(R^ 5 )_fc+i(l) 

and Xk,Xk+i lie on the same curve of W(^R^s),k(t)- Then X^. is the start solution 
iZ{R,s),T, which shows that Xq G I {R, S) G Sols}. 

First set i?(0),... ,R{ro) to be the unique chain from [1,... ,p] to [1,... ,p —l,p-|- 
ro], and similarly for S'(O),... ,S(rQ). Then Xq G hF(_R, 5 ),o(l) and hence lies on a 
unique curve in hF(A;, 5 ),o(^)- Let Xi be the point on that curve with t = 0. By 
Proposition (4), 

2^1 e F{rj(ro),R(ro),Li)( 0 ) = ^ L{R(ro),/l,Ul) ( 1 )' 

PGR{ro)*l 

Let i?(ro-|-l) be the index /3 such that Xi G F(_R(ro),/ 3 ,Li)(l)- Define S'(rQ-l-l) similarly. 
Then Xi G hF(A;^ 5 ')^i(l). 


NUMERICAL SCHUBERT CALCULUS 


19 


In general, suppose that we have constructed -R(O),... , Rij-Q + k), 5(0),... , 5(rQ + 
fc), and Xfc G Then lies on a unique curve in ■ Let 

be the point on that curve at t = 0. Let c be minimal subject to /c < ri + ■ ■ ■ + Tc 
and set a = R{rQ + ■ ■ ■ + Tc-i) and f3 = R{rQ + k). If fc + 1 < ri + ■ ■ ■ + Tc, then by 
Proposition (3), there is a unique index 7 G * 1 such that G If 

k + 1 = ri + ■ ■ ■ + Lc then by Proposition |3.4| (4), there is a unique index 7 G /?* 1 such 
that Xfc+i G Set R{rQ + /c + 1) =7 and likewise dehne 5(rQ + fc + 1). Continuing 
in this fashion, we construct the chains R and 5, and Xj for 0 < j < r. 

We show that R increases everywhere, except possibly at 1, ro + 1,... , tq + • ■ • + 
ra_i+l, and hence R & T{rQ,... ,ra). Similar arguments show that S G T(rQ,... 
which will complete the proof. Suppose/c +1 ^ {1, ri +1,... ,ri + -• • + ra_i + l}. Let 
c be minimal subject to A; < ri + ■ • ■ + rc and let a = i?(ro + • • • + rc_i), /3 = R(ro + k), 
and 7 = R(ro + k + 1). Then by Proposition (3) and (4), /3 7 - The condition 

j{P, 7 ) > j(u, P) in fhe definition of /3 7 ensures that R increases at A: + 1. □ 

4. Homotopy continuation of overdetermined systems 

Numerical homotopy continuation is a method for hnding the isolated solutions of 
a system 

(22) F{X) = 0 

where F = (/i,... , /„) are polynomials in the variables X = (xi,... , x^v). First, a 
homotopy H{X, t) is found with the following properties: 

1. H{X,1) = F{X). 

2. The isolated solutions of H{X, 0) = 0 are known. 

3. The system H{X,t) = 0 dehnes hnitely many (rational) curves crj(A), and each 
isolated solution of (^) is connected to an isolated solution (Ji(0) of H{X, 0) = 0 
by one of these curves. 

Given such a homotopy, numerical path continuation is used to trace these curves 
from solutions of H{X,0) = 0 to solutions of the original system (^^. When there 
are fewer solutions to F{X) = 0 than to H{X,0) = 0, some curves will diverge or 
become singular as A — 1 , and it is expensive to trace such a curve. 

When N = n, the system (^^ is square and the homotopy 

(23) H{XR) := tF{X) + (l-A)G'(X), 

where G{X) = — Oi,... — ajy) with di := deg(/j) and Oj 7 ^ 0, gives 

curves. This is the Bezout bound for a generic dense system F. 

In practice, F{X) = 0 may have fewer than Y\.di solutions and we desire a ho¬ 
motopy with no divergent curves. Methods for such deficient systems which reduce 
the number of divergent curves are developed in [^, |^. When the polynomials 

/i,... , /„ have special forms B 0 . then such homotopies (p3|) are constructed 
where G{X) shares this special form. When the polynomials /i,... , fn are sparse. 
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polyhedral methods give a homotopy. The SAGBI homotopy algorithm (§2.3) 

is in the same spirit. We exploit a special feature of the coordinate ring of the Grass- 
mannian to obtain a homotopy between the system (2.2) we wish to solve and one 
(2.12) whose solution may be obtained using polyhedral methods. Moreover, there 
are generically no divergent curves to be followed. 

The overdetermined situation of n > iV is more delicate. One difficulty is hnding 
a homotopy H{X, t) for an overdetermined system as generic perturbations of F 
have no solutions. In §2], this difficulty is avoided as follows: The system F = 
(/i,... , fn) is replaced by N random linear combinations of the /i,... , /n yielding 
a square system whose isolated solutions include all isolated solutions of F{X) = 0, 
but typically many more. They then find all isolated solutions of this random square 
subsystem. 

For the Grobner and Fieri homotopy algorithms, we gave (in §§2.2 and 3.2-3) 
homotopies H{X,t) = {hi{X,t ),... , hn{X,t)) and solutions <7^(0) at t = 0 as above. 
For these, there are generically no divergent curves. To efficiently follow the curves 
<Ji(t), we select a square subsystem MH{X,t) of H{X,t) (M is an (iV x n)-matrix). 
If the Jacobian of MH at each (7^(0) has the same rank (iV) as does the Jacobian of 
H{X,t), then the curves ai{t) remain components of the algebraic set dehned by the 
equations 

MH{X,t) = 0. 

Moreover, other components of this set meet the curves cTj(t) in at most finitely many 
points f in C — {0}. Thus, we may use the square subsystem MH{X,t) to trace the 
curves (Ji{t) along some path from 0 to 1 in the complex plane. We remark that in 
practice, M may be chosen at random. 

5. Applications 

The algorithms of Sections 2 and 3 are useful for studying both the pole assignment 
problem in systems theory and real enumerative geometry 

We describe the connection to the control of linear systems following [Q . Suppose 
we have a system (for example, a mechanical linkage) with inputs u G R"* and 
outputs y G for which there are internal states x G R"" such that the evolution of 
the system is governed by the first order linear differential equation 

, , X = Ax + Bu, 

^ ^ y = Cx. 

If the input is controlled by constant output feedback, u = Fy, then we obtain 

X = {A + BFC)x. 

The natural frequencies of the controlled system are the roots si,... , Sn of 
(25) (p(s) := - A- BFC). 
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The pole assignment problem asks, given a system and a polynomial <p(s) of 
degree n, which feedback laws F satisfy (pS])? 

A standard transformation (c/. |^, § 2 ]) transforms the inpnt data A,B,C into 
matrices N{s), D{s) of polynomials with det(Zi)(s)) = det(s/ —A) and N{s)D{s)~^ = 
C{sl — A)~^B snch that 


(26) 


(p(s) 


det 


F D{s) 
I N{s) 


Here I is the p x p-identity matrix and the feedback law F is an m x p-matrix. If we 
let 

then F gives local coordinates on Grass{p, m + p) and ( P 6 | ) is eqnivalent to 

Xr\K{si) 7 ^ { 0 } for i = 


and K{s) : = 


D{s) 

JV(s) 


These conditions are independent for generic A, B, C and distinct Sj, hence n < mp 
is necessary for there to be any feedback laws F. The critical case of n = mp is an 
instance of the sitnation in § 2 . 

In P] homotopy continuation was used to solve a specihc feedback problem when 
(m,p) = (3,2). From this result, they deduced that the pole assignment problem is 
not in general solvable by radicals. Despite this success, only few non-trivial examples 
have been computed in the control theory literature . 


An important question is whether a given system may be controlled by real output 
feedback |^, |^. That is, if all roots of ip{s) are real, are there real feedback 
laws F satisfying (p 6 D? Real enumerative geometry asks a similar question: Are 
there real linear subspaces Ki,... , Kn in general position with dimiFj = m + 1 — ki 
and ki + ■ ■ ■ + kn = mp such that all p-planes meeting each Ki nontrivially are real? 


When either m or p is 2 |^, n < 5 |^, or when m = p = 3 and the k^ = 1 |^, the 
answer is yes. In fact, the Fieri homotopies arose from these investigations. 

B. Shapiro and M. Shapiro give a precise conjecture relating both applications. 
Suppose 

Ki{s) := [7(s),7'(s),7"(s),... ,7(”"+^“^*)(s)], 
where 7 (s) is a parameterization of a rational normal (non-degenerate) curve in 
of degree m -|- p — 1. One such choice is 


(27) 


7 (s) = transpose[l, s, s'^ 


m+p-ll 


Geometrically, Ki{s) is the (m -|- 1 — A: 4 )-plane which osculates the curve 7 (s) at s. 
Such osculating m-planes have been used to prove non-degeneracy results in control 
theory. 


Conjecture 5.1 (B. Shapiro and M. Shapiro). Let si,... ,Sn be distinct real num¬ 
bers and suppose Ki{si) osculates 7 at Si and ki kn = mp. Then each of the 






























22 


BIRKETT HUBER, FRANK SOTTILE AND BERND STURMFELS 


finitely many p-planes X C 0™+^ satisfying X fl Kfisfi 7^ {0} for i = 1,... ,n zs 
defined over the reals. 


6. Computational results 


Our algorithms have been tested successfully in MATLAB, hnding all 14 solutions 
in the case {m,,p) = (4, 2) for both the SABGI and Grobner homotopy algorithm, and 
all 15 solutions when {m,p) = (6,2) and ki = ■■■ = kQ = 2 for the Fieri homotopy 
algorithm. 

At present, the SAGBI and Grobner homotopy algorithms have been fully imple¬ 
mented. Some timings from trial runs of these algorithms on a Sparc 20 are displayed 
in Table |l[ The input for these were mp random complex [m + p) x m-matrices. 

We provide a comparison to methods based upon Grobner bases. Table |I| also gives 
the time on the Sparc 20 for the system Singular [|T^] to compute a degree reverse 


lexicographic Grobner basis for the polynomial systems: 


det [X I Ki] =0, for z = 3,... , Kmp- 

Here X is expressed in local coordinates for f2[i3] fl and the K ^,... , K^p are 
[m + p) X m-matrices with random integral entries between —4 and 4. A degree 
reverse lexicographic Grobner basis is the input for some alternative numerical poly¬ 
nomial systems solvers (e.g. eigenvalue methods [§]). We note that the Grobner basis 
calculation did not terminate within one week in the case (m,p) = (6,2). 


m 

p 

d{m,p) 

SABGI homotopy 

Grobner homotopy 

Grobner basis 

3 

2 

5 

<1 

< 0.5 

< 0.5 

4 

2 

14 

47 

6 

19 

5 

2 

42 

373 

408 

149,897 

6 

2 

132 

3,364 

8,626 

CXO 


Table 1. Time (in seconds) 


The hnal version of this paper will include data from implementations of the Fieri 
homotopy algorithm, as well. 
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